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Abstract 

We present tight-binding calculations of the spin torque in non-collinear magnetic tunnel junc- 
tions based on the non-equilibrium Green functions approach. We have calculated the spin torque 
via the effective local magnetic moment approach and the divergence of the spin current. We 
show that both methods are equivalent, i.e. the absorption of the spin current at the interface 
is equivalent to the exchange interaction between the electron spins and the local magnetization. 
The transverse components of the spin torque parallel and perpendicular to the interface oscillate 
with different phase and decay in the ferromagnetic layer (FM) as a function of the distance from 
the interface. The period of oscillations is inversely proportional to the difference between the 
Fermi-momentum of the majority and minority electrons. The phase difference between the two 
transverse components of the spin torque is due to the precession of the electron spins around 
the exchange field in the FM layer. In absence of applied bias and for a relatively thin barrier 
the perpendicular component of the spin torque to the interface is non-zero due to the exchange 
coupling between the FM layers across the barrier. 
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The effect of current-induced switching on the orientation of magnetic moments in mag- 
netic heterostructures has recently attracted significant interest due to its potential applica- 
tions to spin electronics , such as current-controlled magnetic memory elements. When 
a current of spin polarized electrons enters a ferromagnet, there is a transfer of angular mo- 
mentum between the propagating electrons and the background magnetization of the layer. 
The concept of switching the orientation of a magnetic layer of a multilayered structure 
by the current perpendicular to the layers was originally proposed by Slonczewski Q| and 
Berger j^, and has been followed by others 0, 0]. This effect, often referred to as "spin 
torque" , may, at sufficiently high current densities, alter the magnetization state. An alter- 
native mechanism of current- induced switching was put forth by Heide et al 0] and Zhang et 
a/ jsl in which the current across the magnetically inhomogeneous multilayer produces spin 
accumulation which establishes an energy preference for a parallel or antiparallel alignment 
of the moments of the magnetic layers. 

The purpose of this work is to employ a one- dimensional single-band tight-binding model 
and calculate the spin torque in magnetic tunnel junctions consisting of two ferromagnetic 
layers (FM) with non-collinear orientation of the magnetization, separated by an insulating 
(I) thin layer. The calculations are carried out using the non-equilibrium Green's function 
technique in the framework of the Keldysh formalism j^. We have calculated the spin torque 
using two different approaches: the first uses the local magnetic moment of the conduction 
electrons and the exchange splitting of the (i-band; the second uses the divergence of the 
transverse component of the spin current. The results demonstrate that these two approaches 
are equivalent, as hypothesized originally by Edwards et al 

The method for calculating the non-equilibrium Green functions for the non-collinear 
system is a generalization of that introduced by Caroli et al for nonmagnetic leads, where 
one replaces the Green functions with 2x2 matrices in spin space. For this purpose we 
divide the Hamiltonian into two parts: one which describes the spin-average medium and 
the other proportional to the spin splitting of the band structure, i.e.. 

Here, the subscript indexes indicate the atomic sites, the superscript indexes refer to spins, 
and e^^^-* and t^'-^^ denote the spin-dependent on-site energy and nearest-neighbor hopping 
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matrix elements, respectively. The one-electron Slirodinger equation becomes 

/ 



E 

pi 



(^ESpp^ Hpp^I 6Hpp^ 



iJpiq iJpiq 



(2) 



C0S7 siri'^ 
y sm7 — C0S7 

where Qpq"'^ is the retarded Greens function for the isolated semi-infinite leads, 7 is the angle 
between the magnetizations of the FM leads in the plane perpendicular to the current, and / 
is 2x2 unit matrix. Following Caroli ^, we next calculated the retarded Green function for 
the entire system. Finally, by solving the kinetic equation we evaluated the non-equilibrium 
Green function G^. 

Having calculated the G^, the effective local magnetic moment on site i in the right FM 
lead can be obtained from it by 
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where cr is the Pauli matrix vector. The spin current is given by 
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where T""^'""^ is the hopping matrix element in the right lead. 

We model the FM/I/FM junction with = 0.794t^ = 1.79At^, = = 0.4 eV for 
the FM leads, and e = 3th, tb = 1 eV for the insulating layer llj. The insulator consists of 
five layers. We assume that the potential drop is confined to the barrier and is linear within 
the barrier. We have calculated the torque from 



T = Ax 



(5) 



where A = z{e^ — e^)/2fiB is the exchange field along the z-direction, and, alternatively, 
from the divergence of the spin current 



V Q. 



(6) 



In Fig.l we show the spatial distribution of the components of the spin torque on the 
right FM lead that are perpendicular (Tj_) and parallel (T||) to the plane of the layers for 
7 = 7r/2. These are calculated both via the local magnetic moment approach using equation 
(5) (solid and dashed curves) and via the divergence of the spin current from equation (6) 
(solid and open symbols). The direction of eletron flow is from the right to the left FM lead. 
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FIG. 1: Site-dependence of the perpendicular (Tj_) and parallel (Tj|) components of the spin torque 
on the right FM lead calculated via the local magnetic moment approach (Eq. (6))and the diver- 
gence of the spin current (Eq. (7)), for 7 = 7r/2 and ^ = 0.1 V 



One can see that both methods are equivalent giving identical results for the spin torque, i. e. 
the absorption of the spin current at the right FM electrode is equivalent to the exchange 
interaction between the electron spins and the local magnetization. This equivalence is valid 
for an arbitrary angle 7 between the magnetization of FM layers. The angular dependence 
of the spin torque is proportional to sin 7. 

Both T± and T\\ components of the spin torque oscillate with different phase and decay 
with distance in the FM layer from the right FM electrode. The period of oscillations 
is inversely proportional to the difference between the Fermi-momentum of spin-up and 
spin-down electrons. The phase difference between the two components of the spin torque 
is due to the precession of the electron spins around the exchange field in the FM layer. 
The components fi\\ and fi± of the effective magnetic moment are referred to by previous 
studies^, 0] as "spin accumulation". We find that in general, T± and T\\ are comparable in 
magnitude and they can be of the same or opposite sign depending on the position in the 
right FM electrode. 

In Fig. 2 we show the spatial distribution of the parallel and perpendicular components of 
the spin torque at zero and finite bias. At zero bias, Tji = while T± ^ 0, indicating that T±_ 



plays the role of the exchange coupling between the FM leads across the barrier 



121. At finite 



bias, the role of current-induced spin torque is played by T|| and T^iV 0) — T^iV = 0) . 
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FIG. 2: Spatial distribution of the parallel and perpendicular components of the spin torque at 
zero and 0.1 V bias for 7 = 7r/2. 

The exchange couphng explains irregular oscillations of T± in Fig. 1 and Fig. 2 at finite 
bias close to the interface. Since it decays rapidly with the distance from the interface (see 
Fig. 2, T± aX V = 0) , the oscillatory behaviour of T± becomes normal again in the bulk of 
the FM layer. 

So far we have discussed the spin torque behaviour for the case of positive applied bias, 
i.e. when electrons flow from the right FM electrode to the left one. We consider the spin 
torque in the right FM electrode, where the magnetization orientation was chozen to be along 
the z axis. In this case, only those electrons reflected from the left interface are polarized 
transversly to the magnetization of the right FM layer thereby contributing to the current 
induced spin torque. The situation changes drastically when one reverses the current flow, 
i.e. at negative bias. The majority of electrons transmitted from the left FM layer into the 
right one are transversely polarized and lead to a signiflcant increase of the current induced 
spin torque. Both situations are demonstrated in Fig. 3, where Tj| and T± are plotted for 
positive and negative applied bias. 
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FIG. 3: Spatial distribution of the parallel and perpendicular components of the spin torque for 
positive {V = 0.1 V) and negative bias {V = —0.1 V). 
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